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Abstract The sets of the synchronous equations are derived from the sets of non- 
PLh \ synchronous equations The analytical solutions are given by solving the set of differential 

. equations. The results of the evolutionary tendency of the orbit-spin are that the semi- 

major axis shrinks gradually with time: the orbital eccentricity dereacses gradually with 
^ ^ time until the orbital circularization; the orbital period shortens gradually with time and 

I ■ the rotational angular velocity of primary component speed up with time gradually before 

O I the orbit-rotation achieved the circularization The theoretical results are applied to evolu- 

4— > ■ tion of the orbit and spin of synchronous binary stars Algol A, B on the main sequence 

^ I phase The circularization time and life time (age) and the evolutional numerical solutions 

of orbit and spin when circularization time are estimeted for Algol A, B. The results are 
discussed and concluded. 

> 
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1 INTRODUCTION 

. The tidal friction plays an important role in t he evo l ution of the o rbit and sp i n of close binary s ystem. 

^ ■ Earliest, the author researching this topic is IZahnI (1 19651 Il966allbll3. Il975l) . lAlexandeJ (1 19731) firstly 

studied th e dynam i cal pro blem of th e tidal friction in close binary system by using the method em- 
ployed by iDarwi ^ (Il879h . Later on, iHutI (Il980lll98lh generalized the method given by lAlexandeil 
(Il973b . He studied the stability of tidal equilibrium and tidal evolution in close binary system by us- 

■ ing the method of energy and angular momentum. But their research dealt with a few synch roniza- 
\ tion. The sequential research for the synchronization of rotation are given bv IZahnI (Il977l Il978h 



5^ , iRaiamohan & VenkatakrishnanI (Il98ll) ever studied synchroni zation in binary stars . iGiuricin et al.l 

(Il984ah researched synchronization in eclipsing binary sta rs and|Giuricin et alj (Il984bl) also researched 



synchronization in early-type spectroscopic binary stars, IZahn & Bouchetl ( 19891) studied mainly the 



orbital circular ization of la t e-type binar y stars on the pre-main sequence and the theoretical results are 
given based on lZahnl(ll989b . |Panl d 19961) calculated the ci rcularization time scale by using the two mech- 
anisms: one is the equilibrium tidal mechanism given bv'Zahn' d 19771) . another is purely hydrodynamic 
mechanissm given by Tassoul (1987). Keppens et al. (2000) studie d the rotational evolution of binary 
stars system: synchronization and circularization. iHuang & Zen g (2000 ) also research ed evolution of 



non-synchronized binary stars with 9 solar mass and 6 solar mass. Meibo nTet al.l (l2005h studied obsera 



tion tidal synchronization in detached solar-type binary stars and Meibo m et al. (12006 ) also researched 
an observational study of tid al synchronization in solar-type binary starss in open clusters M35 and 
M34. Although the author lLil (1 1 998[ 12004 l2009l) studied some methods forjudging the synchronization 
of rotation of binary stars, but he does not studied the evolution of orbit-rotation of synchronous binary 
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stars. In the present paper the author examined the evolutional tendency of orbit and spin of synchronous 
binary stars on the main sequence phase. 

2 THE EVOLUTIONAL EQUATIONS OF SYNCHRONOUS BINARY STARS DUE TO THE 
TIDAL FRICTION ON THE MAIN SEQUENCE 

The secular evolutional equations of the semi-major axis, a, eccentricity, e, a nd the rotatio nal angular 
velocity, CI, due to tidal frication for non-synchronous binary stars are given bv lZahnl(ll989l) . 
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Where AI and R denote the mass and the radius of primary star, q = ^1' denotes the mass of the 
secondary star, uj denotes the orbital angular velo city (mean motion), I d enotes momentum of inertia, 
the convective friction time and A " are given bv dZahn & Bouchellll989h 
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Here L denotes the luminosity of the primary star. 

Next, one derives the evolutional equations of synchronous binary stars. IZahn & Bouchej (Il989l) 
point out that when the two components rotate with the orbital motion in synchronism: \luj — mi}\ = uj, 
then, all tidal coefficients are identical A''" = A except for A^^. Hence in Equations ([T]) - ([3]) A^^ = 
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A, A ^ X. When we consider the two components rotate in synchronism, Q 



or ^ ~ 1. Substituting these conditions into the Equations ([T]i - Q, the secular Equations ([Til - (O are 
reduced to the following simplified secular synchronized equations 
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We may also write down a supplementary secular equations according to keplerian third law 
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(8) 
(9) 

Where Porb denotes the orbital period and PRot denotes the rotational period 

Substituting the Equation ^ for de/dt into the following equation, we get the circurization scale 

time 

e tf (a' 



^cir- - 21g(l + g)AUJ ■ ^^^^ 

In the following one uses the analytical method to solve the evolutional Equations (|4]l-(|9]l with the 
eccentricity e as an independent variable 

This paper considers the evolutional tendency of the orbit-spin of the synchronous binaries before 
the orbital circularization on the phase of the main sequence. We may assume that the radius, P, of the 
primary star may be regarded as no variation, i.e, P is a constant on the phase of the main sequence star, 
but their separation or semi-major axis is variable due to the action of the tidal friction. 

Combining the Equation Q with Q, we obtain the differential equation 

ij = fe. (H) 

a de 7 

Integrating this equation, we get 

19 

a = aoexp[y(e^ - eg)]. (12) 
Substituting the Equation ( fT2] i into the Equation (|5]l, we obtain the equation; 

= -219(1 + cxp[-i|?(e^ - eg)]. (13) 

e at tf ao 7 

By letting c — 152/7, the differential Equation ( fT3] l can be written 

-expc - el)de = exp -ceg — ^-^ '-de = -21q(l + q) — (—) dt. 

e e tf ^ao 

Using the expansion of series 

exp(ce^) = 1 + ce^ + ^c^e'' + ^c^e*^ + 

Integrating the differential equation above, yields 

exp(-ceg)[hi(e) + ^ce^ + ^c^e^ + = -21^(1 + 9)f - to). 

z o If ao 



We obtain time scale in terms of eccentricity, e, by neglecting the term with 
Here 

A . P. 8 /152 2\ 
g = -(-) exp(— eg). (15) 

Combining the Equation (|4| with we get the equation 

57{l + q),I.dn 
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Integrating this equation, we obtain 

r! = r!o exp [ - — {-^) {—) (a^ - al) 
The from Equation ( fT2] i a ~ ao exp [-^^(e^ ~ eo)] ' s° '^^ — ^xp [2 x ^(e^ ~ e§)] , hence 

2 2 2 I r 38 , 2 2\1 
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We obtain the variation of the angular velocity of primary in terms of eccentricity, e 

r 9 

ri = r^o exp < — - ( 
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The integrations of the Equations (|7ll - (|9]l can be obtained 
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Next, one gives the analytical solutions of the secular evolutional equations with time t as an inde- 
pendent variable on the main sequence. 

For small eccentricity, e, such as Algol A, B, e = 0.015, so ^ce^ = 0.0024, the second term of the 
right hand of formula ( fT4l ) may be neglected., and then, we can obtain the eccentricity decreases with 
time from the formula ( fT4] i 

e = eo cxph21g(l + q)Q{t - to)]. 



- eo = e2{exph42(7(l + q)Q{t - to)] - 1} 
Substituting Equation (|22] | or Equation ( l23T l into Equation (fT2b . we get 

f 19 n 

a = ao exp I y eo{exp[-42(7(l + q)Q{t - to)] - 1} 

a^-al = al jexp | y e2{exp[-42<7(l + q)Q{t - to)] -!}}-! 
Substituting the formula (l22l i into the formula ( fTSl l or the formula dZSb into the formula ( fTSI l, we obtain 

r! = 17o exp |- jexp |ye2{exph42g(l + g)Q(t - to)] - 1}| - ijj . (26) 
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The integrations of the Equations (|7]i - (|9]l can be obtained 
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3 THE ORBIT-SPIN EVOLUTION OF THE SYNCHRONOUS BINARIES (ALGOL A, B ) 

The echpsing binary system Algol {f3 Per) at least consists of three components A, B, C. In fact, there is 
also a massive unseen fourth component D (Hopkins 1976). Algol A (primary) is a main sequence star 
(Bg.V). Algol B (secondary) is a subgaint (g Kq). The separation between A and B is very near mutually 
and the s ystem Algol A, B is re garded as a synchronous binaries due to the action of the tidal friction. 
Based on lGiuricin et alj ( Il984al) the mean rotational angular velocity of primary A is ii = 56 km s~\ 
the mean orbital angular velocity (the synchronized velocity), Vk = 55 km s^^, and based on lTanl(ll985[) 
gave Wsini — 55 km s~^, Wsyn = 55 km s^^. So A and B is a near synchronous binaries and based on 
apparent descriptive method forjudging the synchronization of rotation of binary stars (lLill2004 l2009b 
Algol A, B is a near synchronization binaries. Hence this paper chooses Algol A, B as synchronous 
binaries to calculate the orbital circularization and evlution of the orbit and spin before circularization 
on the main sequence. For the data of Algol A, B, we cite the orbital period (Po)orb = 2.8672 day, 
ao = 14.03 Rq, M = 3.7 M^. m' = 0.81 M^, , R = 2.74 R = 3.60 Rr^, q = M'/M = 0.22 
Te = 12010" K dBrancewicz & DworakI [l980l) . eo = 0.015 (iTomkin & Lambert! 1 19781 iHarringtonl 
11984 . no = v/R = 2.9365 x lO''^ rad s'^ a ndtf = tfp.( M/Mc^y/^(Tj^,p^)-'^/^, tfQ = 0.433 yr, 
Te© = 5770" K, tf = 0.2519 yr, L = 2.2 ('Popper"l980). 

Where T p is the effective te mperature (jZahn & Bouchet 19891) 

Recently, Yang et all (|201 ll) present XMM-Newton observation of the eclipsing binary Algol. Their 



results of research ar e valuable. 

As lZahn & Bouchet Cl989h showed that when the coefficient A are all equal, then 

A — k2 (the apsidal motion constant). (30) 

A = fc2 (the apsidal motiom constant) is calculated from the formula given bv lCowling| ( ll938h and 
letting fci = fc2 = fc 

(f)^(l + 16i^) + (f)'^(l + 16(^)- ^ ^ 

Here denote s the period of advance of the apsital line, P' ~ 2.476 year for Algol A, B given by 
lHegediisl(ll988h . Substituting Porb, P', o,, M, M', R and R' into the above formula, we get 

A = fc = 0.003308. (32) 

The moment of inertia / 

/ = KMR'^. 

K is culculated from the formula ~ + f) (ISchatzmanll 1 9631) . The polytropic index n = 3 for 
the main sequence star (Algol A). So if = 4/33 = 0.1212. 

M M 1 



0.02268 X lO^'*" kmr\ exp(152e5/7) = 1.00489. 



/ KMK^ KR^ 
Substituting the values of fc, tf,R,a and exp(152eQ/7) into the formula ( fTSl l, we obtain 

g = 2.7 X 10-^yr-^ (33) 

Let us estimate the numerical solutions when the orbit of Algol A, B achieves the circularizatuion 

At first, we evaluate the circularization time scale Substituting the values of q, Rq, oq, tf and 
A = fc = 0.003308 into the formula ( fTOl ). we obtain the circularization time scale 

t^i, 6.4589 X 10*^ yr. (34) 

Next we estimate the numerical solution of evolutionary tendency of the orbit-spin when Algol 
A, B achieved the orbital circularization time. By letting the intial time to = 0, and substituting 
t^i, = 6.4589 yr into the formulas ^ or ^16^-^, we get a = 14.O226i?0, e = 0, 0056, 
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w = 2.1931 rad d-\ Porb = 2.8649 d, ^2 = 2.5536 rad d-\ PRot = 2.4605 d; 5a = -0.0074i?o, 
(5e = -0.0094, Su = +0.0017 rad d-\ ^Porb = -0.0023 d, 5n = +0.0156 rad d-\ SP^ot = 
-0.0160 d. 

The life time (age) is based on stellar mass-lose M, i.e 



M 

tnic = — • (35) 
M 

The value of its may be calculated from the formule given byi owers &Deemingl(ll984 

or calculated from the formula given bv lNieuwenhuiizen &de Jaged(ll990h 

log M = -14.02 + 1.24 log(i/Le) + 0.81 \og{R/Re) + 0.16 log(M/Af0) (37) 
Substituting these datd of M, R and L into the formula (l36T l and (l35T l, we get the life time (age) 

tiif^ = 7.5703 X 10^ yr (age). (38) 
The lime time for the speed up of spin 

tn = ^= tf/36q^Ake\\-)^ = 4.2447 x 10« yr. (39) 

The orbital spiral time (the collasped time of system) 

ta = ^= tf/lUq{l + q)ke^{R/af = 5.2273 x 10^ yr. (40) 

4 DICUSSION AND CONCLUSIONS 

(1) The set of the Equations (|4li-(|6l) fits in synchronous binaries on the pre-main sequence, main se- 
quence and post-main sequence phase according to the radii of a stars. The raduis of the late type 
star is variable due to the gravitational contraction on the pre-main sequence phase. The raduis of 
gaint star is variable possibly due to the expansion of shell on the post-main sequence phase. On the 
main sequence phase the raduis of star is stable. Its raduis can be regarded as constant. These refer 
to the raduis of primary star because in the Equations dU-® R denotes the raduis of the primary 
star It does not refer to the gaint star 

(2) The research of this paper differs from that of IZahn & Bouchetl (Il989l) in some aspects. Zahn & 
Boucher's paper researched the orbital evolution and circularization of binary stars on pre-main 
sequence phase by analytical method and for non-synchronized equations by the method of nu- 
merical integration. In the analytical method the radii of binaries are variable due to gravitational 
contraction, but the semi-major axis is not variable including the main sequence phase. However 
the present paper studies the orbit-spin of binary stars on the main sequence phase by analytical 
method in which the star's raduis is not variable and the semi-major axis is variable due to tidal 
friction. In Zahn & Boucher's paper they must use the method of numerical integration to solve 
non-synchronized equations. However in the present paper the author may use the analytical method 
to solve the synchronized equations. Zahn & Boucher estimate that the eccentricity decreases from 
0.005 to 0.0043 in 10 billion for binaries 0.5 Mq + 0.5 Mq on the main sequence phase. The 
present paper estimates that the eccentricity decreases from 0.015 to 0.0056 in 6.45 mega for bina- 
ries 3.07 Mq + 0.81 Mq on the main sequence phase. Hence both results are different. 

(3) The results of the solution for integrating diferential equations by analytical mathod are some less 
different with that by the method of numerical integration. For example, the semi-major axis a = 
14.O226i?0 for Algol A, B when the circularization time (6.4589 x 10^ yr) by the former method 
and a = 14.O126i?0 by the latter method. However this difference is very small. 
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(4) The circularization time occurs before the life time (age), and the lime time for the speed up of spin 
and the spiral time (the collasped time of the system) occur after the life time (age). Hence the latter 
both are meaningless. 

(5) In the system of Algol A, B, C the tidal frictio n in triple stars jKiseleva et al.|[T998l) and the perturb- 
ing effect of the third star (Algol C) ( lLill2006b may decircularize the orbit of the second star (Algol 
B). In this paper one does not consider these effect. 

One obtained the following conclusions: 

(1) The eccentricity decreases gradually with time until the orbital circularization, i.e until e decreases 
to the circularization time. 

(2) The semi-major axis shrinks gradually with time or with eccentricity decreases. 

(3) The orbital period shortens gradually with time or with circularization. 

(4) The rotational angular velocity of primary component speeds up with time gradually. 
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